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22 Abstract 

/-^ We argue that in cold, dense quark matter, in the limit of a large number of colors the ground 

^^ state is unstable with respect to creation of a complicated Quarkyonic Chiral Spiral (QCS) state. 



o 



X 



in which both chiral and translational symmetries are spontaneously broken. The entire Fermi 
surface is covered with patches of QCSs, whose number increases as the quark density does. The 
low energy excitations are gapless, given by Wess-Zumino-Novikov-Witten model plus transverse 
kinetic terms localized about different patches. 
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I. INTRODUCTION 

In cold, dense quark matter, one expects several phenomenon: a transition to a deconfined 
phase, color superconductivity [1], and for a large number of colors, N^, chiral density waves 
[2-5]. 

Recent theoretical work has argued that at least for a large number of colors, A^c, further 
phenomena can occur [6-12]. There can be a phase in which the baryon density is large, 
but that is nevertheless confining. While such a "Quarkyonic" phase is certainly present at 
large A*",,, it also exhibits several properties which may have been observed in experiments 
at the CERN SPS [9]. 

Ref. [10] suggests that in a Quarkyonic phase, the local breaking of chiral symmetry 
is driven by the appearance of chiral density waves, which were termed Quarkyonic Chi- 
ral Spirals (QCS). This arises from the condensation of particle-hole pairs near the Fermi 
surface, and is similar to the mechanism of chiral condensation in the vacuum. The essen- 
tial difference is that in the presence of a Fermi sea, the dominant condensation involves a 
particle-hole pairs with nearly opposite Fermi momentum. The net momentum of the pair is 
then nonzero and large, so that not only chiral symmetry, but also translational symmetry, 
is spontaneously broken through the formation of a chiral condensate which was termed a 
(Quarkyonic) Chiral Spiral. This condensate rotates between the usual chiral condensate, 
ipip, and a helicity condensate, ip'^o'fz'^^- This is closely analogous to a chiral spiral in 1 + 1 
dimensions, which rotates between ipip and ipj^ip [13]. 

The analysis in Ref. [10] was restricted to chiral spirals which form in one direction. In 
this paper we argue that the entire Fermi surface is covered with patches of QCSs, where 
each patch corresponds to a different direction for a single QCS. The long-range forces 
present in Quarkyonic matter imply that the gaps are strongly momentum dependent. As a 
consequence, at large angles the interaction between particles and QCSs is suppressed, and 
QCSs at large angles to one another interact weakly. 

We estimate the width of the patches, and find that the number of patches increases with 
the density. As the density increases, the Fermi surface is characterized by increasingly high 
symmetry, and a series of phase transitions occur. All of these transitions occur in a confined 
phase with locally broken chiral symmetry, and thus are distinct from the transitions for the 
restoration of chiral symmetry, or for deconfinement. 



Ref. [10] showed that QCSs dominate at intermediate chemical potentials, up to ;U ~ 
80 GeV. At higher fi, the effects of perturbative interactions must also be included, and also 
generate chiral spirals. It was previously suggested in Refs. [3] and [4] that perturbative 
chiral spirals also cover the Fermi sphere in patches. Our analysis extends and compliments 
theirs. 

II. THE QUALITATIVE PICTURE FOR THE CONFINING MODEL 

In the limit of large A'^c, the gluon propagator is unaffected by the quarks. This means 
that at low temperature, the gluon propagator is the same as in the confined vacuum. We 
thus use the effective gluon propagator 



-^, Dik{uj,q) = — — - 



^oo(c^,q) = ^^, Afc(w,q) = 2 I 2 ■ (1) 



This propagator is valid in Coulomb gauge, diAi = 0, at small momenta, \q\ < Aqcd- This 
propagator is very long range, and through the exchange of Aq fields, strongly enhances the 
density-density interaction at small momentum transfer. This reflects the linear confinement 
of mesonic excitations. Following [10], we concentrate on nonperturbative contributions from 
Dqq and ignore perturbative effects from Dik. 

We summarize the results of [10], and give a qualitative interpretation which we use to 
construct a Fermi surface with multiple patches. Taking the gluon propagator Dqq of Eq. 
(1), the Schwinger-Dyson equation for the quark self energy can be constructed, especially 
for excitations near the Fermi surface. It is 

X:(p) + S„(p) = - / ^ ilotA) S{k; S) (7otB) /^o1f (p - k) . (2) 

Near the Fermi surface, pz ~ fi, and \pt\ ~ 0. Since the dominant contributions come from 
interactions with small momentum transfer, the integral is dominated by virtual quarks 
within ~ Aqcd of the edge of the Fermi sea. In this region the transverse momenta of the 
quarks can be ignored, and the integral equation factorizes. 



dk.QdkzCp'kT S{kQ,kz,kT) Dqq[p — k) — )■ dkodk^ S{ko,kz,OT) dkx Dqq{p — k) . (3) 



Ot)J' 



Integrating over the gluon propagator over kx gives the Schwinger-Dyson equation for QCD 
in 1 + 1 dimensions in axial gauge, A^ = 0. 



When spatial dimensions are compactified, usually motion in the small, compact dimen- 
sions is energetically disfavored. In contrast, here the quarks move easily in the transverse 
dimension, with essentially no energetic penalty. Thus counting the number of states with 
different transverse momenta matters, and is taken into account in the second integral of 
Eq. (3). This smears the four dimensional, confining gluon propagators, and reduces them 
to propagators in two dimensions. 

Next consider the particle-hole correlations. These can be treated in a Bethe-Salpeter 
equation for quarks and quark holes, matching north and south poles of the Fermi sea. 
By the same reasoning as applies for the Schwinger-Dyson equation, the Bethe-Salpeter 
equation also reduces to that for QCD in two dimensions. 

One might ask why we only choose states at opposite ends of the Fermi sea, and do not 
simultaneously take into account contributions from the entire Fermi surface. For example, 
in the usual analysis of (s-wave) BCS superconductivity, the Coulomb force is additive, and 
one electron feels a sum of interactions from all electrons over the entire Fermi surface. 

This is not true in a Quarkyonic phase. For a system with quark and quark- holes, 
because of confinement what is dominant is the formation of mesons between quarks and 
quark holes. The resulting mesonic states then interact weakly, via residual meson-meson 
interactions ~ 1/A^c- Therefore it is reasonable to first identify the most tightly bound pairs, 
treating interactions between the mesons as a small perturbation. 

This is why we first consider a problem with two patches, and then extend this analysis 
to multiple patches. The most relevant pairing for two patches is the density wave channel 
in which quarks and quark holes move in the same direction with total momentum, ~ 2/i, 
while keeping their relative momentum to ~ Aqcd [10]. Confinement then drives a gap 
~ Aqcd, and is not very sensitive to the quark density. This will allow us to easily estimate 
the transverse size of the patches. 

Finally we briefly mention the mapping between quark bilinears in 3 + 1 dimensions 
to those in 1 + 1 dimensions. The basic point is that the transverse momenta is negligible 
relative to ~ /i, so that quantities which couple to transverse momenta are subleading. Thus 
7t can be dropped, as can terms which mix spin. The spin symmetry converts to a doubled 
flavor symmetry in 1 + 1 dimensions. For a dictionary of quark bilinear operators, see [10]. 



III. THE SINGLE PATCH PROBLEM: QUALITATIVE ARGUMENTS 

In condensed matter physics, usually density waves form because of a nesting of the Fermi 
surface, with particle-hole pairs near the Fermi surface connected by a common wave vector 
[14]. In one spatial dimension, such nesting is automatic, because the boundary of the Fermi 
sea are just two points. This is not true in higher dimensions. For quarks in 3 + 1 dimensions, 
the boundary of the Fermi sea is spherically symmetric, and nesting cannot occur without 
deforming the Fermi sphere. 

To develop insight into how the chiral condensate forms, we first assume that the single 
particle spectrum for quarks is a well defined notion. We then apply mean field theory, and 
ignore all of the subtleties from confinement. In the last part of this section, we argue more 
carefully, albeit qualitatively, as to why this might be correct. 

We first consider an order parameter with only one Fourier component. We use a mean 
field approximation where the chiral spiral acts on quarks like a periodic background, ~ 
Aq cos(Q-r). This problem is mathematically equivalent to that of noninteracting electrons 
propagating in a crystalline periodic potential. The dispersion relations for the particle and 
hole branches are given by 



^^ ^ ^(P)+^(P + Q) ± ^^(P)-^(P + Q)' + |A,|. (4) 

where e(p) is the quark energy. Aq is the magnitude of the condensate induced by the chiral 
spiral waves at a momentum Q, and is treated as a variational parameter. An optimum |Q| 
should be about 2pF, so we can safely replace Aq with M = A2pp. 

Let us establish what the requirements are for Q and the patch size. The 2;-axis is taken 
to penetrate a center of two patches connected by Q, and the momenta p is measured from 
Pf^z- Expanding the energies about the Fermi surface, and taking Q ~ 2pp, p^ ~ Pf, 

p\ 

-r ..., 

e(p + Q) = JiQ-Pz-ppy+pl -PF^Q-2pF-p. + ^ + ... . (5) 

^Pf 

Substituting this into (4) gives 



e(p) = \/{Pz + PfY + pI - Pf ^ Pz + 7^ — + 

ZPF 
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Pf-Pz] +M2. (6) 



The dispersion relation for a single quark can be used to estimate the size of the patches. 
Energetically, the system benefits from opening as large a gap as possible. The edge of a 
patch occurs when a single quark becomes gapless. This refiects the balance between the 
condensation energy and the kinetic energy. 

Thus Q, and the size of the patch, is determined as follows. A lower bound for |Q| arises 
from the condition that the positive branch has nonzero energy, E+ > 0. The minimum 
value of \Q\ occurs when E^ = 0; this happens for p± = 0,pz = Q/2 — pp, which satisfies 

^f° = I - PF + M > 0. (7) 

Thus the ordering vector is not exactly 2pF, but a slightly smaller value, 2{pp — M). The 
minimum value for Q is Qmin = 2{pf — M), caused by the deformation of the Fermi surface 
(see Fig. 1). 

The transverse size of the patches is determined by the condition that the lower branch 
has positive energy. E^ = occurs when 

p^ = 9._pp = _M , ^^^ = M. (8) 

^ 2 ^ ' 2pF ^ ' 

For later use, we refer to p™^^ = ^2pfM = Mj_. 

The main origin of the gap is confinement. In the Quarkyonic phase at large N^, the 
string tension is independent of density, and thus so is the gap. 



M ~ v^ ~ Aqcd • (9) 

Why can we use the quark dispersion relation to estimate the transverse size of the 
patch? After all, with a sharp infrared cutoff for the confining gluon propagator, the quark 
propagator has a pole which diverges as the cutoff vanishes. This divergence can be avoided 
with other infrared regulators, but it is worrisome. 

The point is, however, that any quark excitation is inevitably accompanied by a quark 
hole, with the two forming a mesonic type excitation. For such wavefunctions, infrared 
ambiguities disappear. Consider, for example, the theory in vacuum. All quarks are confined, 
but it is sensible to speak of the constituent mass of a quark, approximately one half the 
mass of a meson. For the mesons near the Fermi surface, then, we can also speak of a 
constituent "mass" . The dependence of this constituent mass on the transverse momentum 
should be given, approximately, by our naive analysis above. 

6 




FIG. 1: The Fermi surface (FS) with three QCSs. Bringing the opposite pieces of FS demonstrates 
imperfect nesting. The spherical parts of the FS on each edge are approximated by Eq.(21). 

Of course, the patch size can be rigorously computed by considering color singlet quanti- 
ties. In 1 + 1 dimensions, this can be done using non-Abelian bosonization. Detailed analysis 
shows that the width of the patch is ~ p™^^, up to a number of order one. 

IV. MULTIPLE PATCHES 

In the previous subsection, we considered one patch problem in (3 + 1) dimensions. The 
size of the patch was estimated by computing the transverse momentum at which the gap 
closed. In reality, the gap remains nonzero, due to interaction with different patches. Since 
the size of the gap is related to that of the order parameter, different patches then interact 
with one another. 

However, a significant simplification takes place when we take into account that the gap 
depends upon the virtuality of the quark. In vacuum, it is known that quarks with high 
virtuality, or in the deep Euclidean region, do not feel the effects of chiral symmetry breaking. 
This suggests that the scattering between a particle and a QCS is small if the angle between 
the two is large. 

To see this more explicitly, note that in the presence of a QCS, the mass self-energy is 

where |Qj| ~ '^Pf- Consider the process whereby a quark, nearly on-shell with momentum 
Po — Pf and p ~ Pf^z, scattering of a QCS with momentum Qi, and thereby acquires 



momentum p — Qi. QCSs are static, so the energy po remains the same. The mass function 
M(p; Qi) contributes if the scattered quark is on-shell within a region ~ Aq(2;j3, 

pI-{p- Q? - -4pMl - COS0) ~ -A|cD, P Q = 2plcose. (11) 

Thus the angle to the 2;-axis, 6, is hmited to \6\ < Aqcd/pf- Physically, this is the same 
condition as for the collinear scattering of two quarks. This condition is special to the long 
range interactions in Quarkyonic matter. 

Hence the mass term in (10) is dominated by patches whose angle is within or tt of the 
incident quark, 

E J (^^(^ - Q^)Mip■, Q^Mp) X e f A^cD - |2p ■Q^-Q-\]+ (h.c). (12) 

In the following, we will take into account interactions between nearest neighbor patches in 
momentum space. 

From Eq.(8), the number of patches is 

where n is the quark density per single flavor. Thus the number of patches increases slowly 
with the density. Since the number of patches is an integer, it must change discontinuously. 
In Fig. 2 we depict a succession of four different coverings of a sphere, corresponding to 
different arrangements of QCSs. At the smallest density there are three mutually perpendic- 
ular QCSs, which form a cube. As the density increases, extra patches emerge at the corners 
of the cube, giving rise to four more QCSs, the blue shape in Fig. 2. The wave vectors of the 
new QCSs are aligned along the diagonals of the cube, so that eight new patches emerge at 
the corners of the cube. Increasing the density further leads to extra patches at the corners 
of the blue shape, leading to the red shape, thence to the green shape, and so on. 

Instead of the above construction, what if one takes smaller patches, say of size Aqcd, 
with a corresponding increase in the number of the patches? This is not possible, though. 
Because the gap arises from a confining force, the size of Q is fixed to be \Q\/2 ~ pp — M. By 
Luttinger's theorem, the total volume of the Fermi sea is conserved, so any volume squeezed 
in the p^-direction must be compensated by an extension in the transverse size, until it 
reaches ~ M_|_. We expect that the presence of other patches slightly reduce the transverse 
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FIG. 2: A succession of QCS configurations corresponding to the different quark density. Black 
lines show a simple cube. 

size, aM_|_ (Aqcd/^± <^ a < 1), while \Q\ slightly increases from the value for one patch. 
Below we will determine a and |Q| by regarding them as variational parameters. 

To match different patches onto one another, consider the part of the Fermi surface where 
several patches meet. From Fig. 2 it follows that the most typical configuration is where 
four patches meet at one corner. We want to establish the conditions for the emergence of 
next patch. Assuming that there are many patches, the curvature at the corners can be 
neglected. We also neglect the influence of all other patches except for those four which are 
adjacent. Then the problem reduces that of a particle subject to the potential generated by 
four QCSs. If Q is a vector in the direction of the corner, then the four wave vectors of the 
patches are Qj = Q ± aM^ej, where 



E' 



0, 



1, Q-ei = 0, M^ 



Q\M <^ \Q\ 



(14) 



with M_|_ = yJ2ppM is the transverse size of the patches, from above. We expect that having 
several patches does not change this, so we take a ~ 1. In high density, we can set Q ~ Q 
with IQI = 2pF, and Q • Cj ~ 0. 
Near the corner, the total gap is 



AQ(r) 



M 



E- 



,i{Q±aM^e,)-r 



Me'' 



cos(aM_|_a;) + cos(aM_Lj/) 



(15) 



4 ^^ 2 

i 

By the previous arguments about virtuality, only particles within ~ Aqcd of a corner con- 
tribute. Away from the corner, where \p\_\ ^ Aqcd, we reduce to the previous problem of 



one patch. 

Since the transverse modulation occurs slowly, we can calculate the spectrum in the adi- 
abatic approximation, separating the motion along Q from the transverse motion. Initially 
we consider the transverse coordinates as constant parameters. Then the energy eigenstates 
yield the two bands of the one patch problem, Eq. (5). Next we consider the energy for the 



effective Hamiltonian, H^rr, describing the transverse motion of the upper (+) and lower 



'—) bands. 



Hf^^i^-PF 




Q 

2 



|A(rx)M = M 



cos(q;M_lx) + cos(aM_L?/) . (16) 



This patch configuration becomes unstable when the gap between the upper and the lower 
bands close. This occurs first when k^ = Q/2. Eq. (16) is the Schrodinger equation 
in a periodic potential. Its solutions and eigenvalues Ei{p±), I = 1,2,..., depend on 
the continuous variable p_L (quasi-momentum) with values confined to the Brillouin zone 
—aM±/2 < {p^,pY) < aM±/2. The solutions are arranged so that Ei < E2 < E^ < .... The 
given QCS configuration becomes unstable if the lowest band, / = 1, crosses from negative 
to positive energy at some value of quasi-momentum. 

Taking into account (14) and making the substitution r_L — )■ {aM±)^^r±, we get the 
Schrodinger equation for the transverse motion: 

(-«VI + aJ)±^^^^±^.±(f-p.))*.(f)*, 

mri + 2™) = exp (^^^^p^) mr) (17) 

where a = {n,m) is a lattice vector. Remember that Q/2 — pi? < 0, so positive and negative 
branch conditions are not automatically satisfied. We have to choose appropriate regions 
for a and |Q|. 

Let us find such a region. The transverse equation can be solved analytically at small 
and large a. At a < 1, we can approximate the potential as a sum of harmonic wells. For 
the negative branch, 

^- ^ ^^ Pf] -2 + 2a-^ + 0{a') 



M M V2 y 4 

-l(cosf^)+cosf^)) expf-^) (18) 

10 



Since {Q/2 — pp) ~ — M, the last term is small. At small a, Eq ^ —3 + 0{a), and satisfies 
the negative branch condition. 

At large a ^ 1, the spectrum E_ is dominated by the kinetic energy, E_ ~ a"^, and is 
positive. Thus the spectrum crosses at some positive value of a ~ 1, in agreement with the 
rough estimates of Eq. (8). When a reaches the critical value, the gap between the bands 
shrinks to zero at some point in the Brillouin zone, creating conditions for the emergence of 
new QCS with a wave vector at this point. 

V. EFFECTIVE ACTION FOR THE COLLECTIVE EXCITATIONS. 

We now derive the effective action for collective excitations. This is possible to do even 
without detailed knowledge of the quark correlation functions. We begin with the collec- 
tive modes located on a particular patch, treating the curvature of the Fermi surface as a 
perturbation. 

Consider two opposite patches of the Fermi surface connected by vector Q, taking the 
2;- axis along Q. Diagonalizing the Dirac Hamiltonian for quarks close to the Fermi surface. 

Ho = f dz(fxM{~R^d,R + L^d,L) + ^{V±R^V±R + V±L^V±L) + ---\ ; (19) 

where R, L fields represent right-moving and left-moving particles along 2;-direction. We 
suppress all indices, including those for spin. We consider A^f flavors in 3 + 1 dimensions, so 
that fermion fields in one-dimensional chain effectively have A^£ = 2N{ flavor indices. 

We treat the system as many one-dimensional chains which couple through transverse 
kinetic terms. We discretize the transverse dimensions, so that the transverse degrees of 
freedom emerge as a result of interchain tunneling [16], 

Ho = J2 [ ^4^( - R\r)dzR{r) + Lt(r)9,L(r)) 

r -^ 

-t^ [i?"^(r)i?(r + a) + L'^(r)L(r + a) |. (20) 

a 

Here a are vectors connecting nearest neighbors, and t and a are chosen so that the spectrum 
of the discrete model approximates the continuous theory. The requirement of nearest 
neighbor tunneling is not necessary, and can be generalized, if need be. Choosing a square 
lattice, the resulting dispersion relation is 

e±{p±) = — 2tN^cos[a • p±] = —2t[cos{pxa) + cos{pya) — 2] (21) 
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In the continuum limit this dispersion relation reproduces p\/2pp on the interval set by (8): 

ta^ = l/2pF, 8t = M. (22) 

Therefore a~^ = {pfMY^'^/2 ~ M±. As we have shown in the previous subsection, the trans- 
verse dispersion can be substantially renormalized by the interaction between the patches, 
Eq. (18). 

The Hamiltonian for each patch is ($ = -R + L) 

H = J2HiDir)+Tt^nn (23) 

r 

Hid = j dz[i{ - R\r)d,R{r) + L\r)dMr)) 

+ f dw^ ($^tA$(r, z)) \z-w\ ($^tA$(r, w)) } (24) 

Ttunn = -tYl [RHr)R{r + a) + L\r)L{r + a)] (25) 

It is important that the interaction is essentially one- dimensional and is concentrated inside 
of each chain. One-dimensional models (24) can be treated by non-Abelian bosonization, as 
was done in [10]. The long range current-current interaction generates a gap of magnitude 
M in the color channel. The remaining gapless modes are described by the Uf^i]^^{l) and 
SUnXN'^) Wess-Zumino-Novikov-Witten (WZNW) models. 

Some care must be taken in the normalization of quark fields, originally in 3 + 1 dimen- 
sions, when we apply bosonization rules to the corresponding fields in 1 + 1 dimensions. A 
one dimensional chain has resolution ~ a^^ ~ M_|_ in the transverse direction, so it is nec- 
essary to multiply dimensionless fields by a~^ when we include the transverse dimensions, 
z — )■ (-2,t). For instance, 

M{z, r) -^ a-'^ X Cxd{U tr[^] tr[/i] + (/i.e.)) (26) 

where the matrix valued fields ?7, (?, and h on the right hand side are boson fields for f/(l), 
flavor, color sector. All of these fields are dimensionless. From the analysis of QCD in 1 + 1 
dimensions, Cxb has dimensions of mass, ~ Aqcd- 

Next we treat Eq. (25) as a perturbation on the WZNW models. Turning on the trans- 
verse interactions, two things happen. While the single particle gap closes at the border 
of each patch, as discussed in the previous section, the gap is restored when we take into 
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account the fact that several patches meet at that point. In the mean field approximation, 
the size of the patch is given by Eq. (8). A more rigorous approach would be to use the 
exact ID quark's Green's function for 1 + 1 dimensional QCD, and to use the Random Phase 
approximation (RPA) expression for the single quark Green's function, in the spirit of Ref. 

[16]: 

GRR = [G-^lM^,P.)-eip^)]-\ (27) 

where Gfi,r,id is the Green's function for right-moving quarks in the one dimensional Hamil- 
tonian of Eq.(24). A similar expression, with R replaced by L, holds for the left-moving 
quarks. Then the size of the patch is determined by the transverse momentum at which the 
Green's function (27) acquires a singularity at u; = 0: 

^ = minG^J,^i^(0,p,) (28) 

At present we do not know the form of Gid', if we assume that 

1 ^(pi 



GMO,P.) = -f[^jf,] (29) 



then we find that 



I3M , /3 ~ 1 . (30) 



2pF 

This agrees with the mean field estimate of Eq. (8). 

The second effect of the interchain tunneling is that it generates a coupling between 
collective modes inside of a given patch. This occurs through virtual scattering processes, 
into states above the single particle gap. In the interaction picture for Hid, time ordered 
product in second order in t generates the interchain coupling, 

t^ J2 f d'^z(fw[R^iz,r)R{z,r + a)][L^{w,r' + a')L{w,r')] 

r,r' ,a,a' 

= -t^a^ ^ / dPzdPw [R]^{z, r)Lf,c' {w, r)] [l\,^, {w, r + a)Rfc{z, r + a)] + ■ ■ ■ , (31) 

r,a 

where /, c express flavor and color indices, respectively. We neglect terms where r ^ r' or 
a ^ a', which do not contribute at this order. 

We integrate over the color sector, leaving color singlet operators. We use coherent fleld 
expressions, L = {a-^cH^) ^l^l/^lc and R = {a~^ClJ^) ^r^r/^rc, which satisfy ^{j^l = U, 

13 



^Rf^Lf = 9ff', and Crc^Lc' = h^c'- Due to f/(l)-flavor-color separation in Hid, we can 
evaluate each sector separately. Note that only color sector has the scale ~ Aqcd ~ M. 
Further, correlation functions in the color sector damp rapidly, while those for the U{1) and 
flavor sectors have power law correlations. This scale hierarchy allows us to use the following 
approximation: 

^ CiD (c] (^\c(.\\ fci 



a? 



C 



ID 



a? 



^ff'i^) {^cc'g{z -W) + {tAJcc'fAiz - W)), (32) 



where Aff/ =: e** : : Qff :. Therefore Eq.(31) becomes 



iN, { ec- 

M2 "" 



( - ^) E / ^'^ ^//'(r' ^)^/'/(r + ^' ^) + om ■ (33) 



The first factor arises from integration of correlation functions in color sector, and 7 ~ 1. 
We have used the fact that correlation functions in the color sector contribute to integrals 
only within a distance ~ M~^. Corrections 0{N^) arise from color non-singlet operators. 
We now have the effective Lagrangian for interchain coupling: 



ia)^ 



J2 f d^ztT[A{r,z)A^{r + ai,z)] = ^^fl^'' f d^x tT[d±A{x)d±A\x)] ,(34) 



where we used 8t = M. For the unperturbed Hamiltonian, the correlation function of A 
scales as 

(A(r,2;)A(0,0))^,, = [^^^^^J . (35) 

Here M serves as the ultraviolet cut-off. The scaling dimensions of the primary fields in the 
WZNW model are [17]: 



0?A 



1 _ iv; - 1/iv; 

Tl ~^ 



(36) 



Therefore the scaling dimension of the operator in Eq. (34) is less than two: that is, 
transverse hopping has generated relevant operators. 

Near zero temperature, the system exhibits long range correlations. The massless modes 
inside of each patch couple strongly to one another, with an energy which depends on 
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their transverse momenta. Their three dimensional dynamics are described by the effective 
Lagrangian 






N[N^PfM 



Pf 



SUiN^) _ KpfM r !/i^,.rn, .n, .tl 






Cwzw[9] + ^tT[d^9d±g^ , (37) 

Pf 



where 



Cwzwig] = ^tr[dLgdLg^ + ^^ /^^ tryd^gWd''g){g^d'g)] . (38) 

Here Ol express derivatives for time and longitudinal directions. We have omitted corrections 
0{M'^/pjp, l/Nc), with the constants rj, rj' ~ 1. 

In one spatial dimension, due to infrared fluctuations there is only quasi long range 
order. After including transverse hopping, our gapless modes become three dimensional, 
Eq. (37), which drastically decreases the role of infrared fluctuations. Consequently, in three 
dimensions there is true long range order, and a well defined ground state. (In contrast to 
Ref. [18], this is only seen when transverse fluctuations are included.) 

Another consequence of the increase in the spatial dimensionality is that nonlinear terms 
in the WZNW action can be treated perturbatively. Goldstone modes of a QCS have a 
spectrum 

a;2 (p) = (n ■ p)2 + K (^j [n x p]^ , « ~ r], r]' . (39) 

The unit vector n is the wave vector of the QCS, with momenta measured from north or 
south pole of the Fermi surface, ppn. This expression is valid for |p_l| ^ M±. It is not 
valid close to the boundary of a patch, where one cannot expand in t, and effects from other 
patches must be included. This shows that at high density the spectrum is less sensitive to 
to the transverse momenta, in accord with Ref. [10]. 

Finally we mention interactions among quantum excitations. Since Eq.(37) has a decay 
constant of 0{Nc), higher order couplings of quantum excitations always accompany sup- 
pression factor of 0(1/A^c), in a similar way as chiral Lagrangian. Therefore, in large A'^c, 
dynamics close to the Fermi surface is described by free Goldstone bosons. The interaction 
between coherent excitations are interesting and will be discussed elsewhere. 
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FIG. 3: The schematic picture of collective excitations for the Fermi surface with three QCSs. 
The magenta ellipses of various intensity correspond to surfaces of constant i^(p). The dark areas 
represent regions of momentum space where the Goldstone modes are subject to decay into quarks. 
The green circle is the original FS. The red arcs represent the portions of FS where spectral gaps 
emerge. 

VI. CONCLUSIONS 

In a Quarkyonic regime, where quarks are confined, instead of treating the entire Fermi 
surface, it is more natural to first construct mesonic degrees of freedom within a particular 
momentum region. Since gaps arise from confinement, the gap energy is almost independent 
of density. 

We analyzed multiple patches of QCSs based on results for a single patch. In going from 
one to multiple patches, it is essential to remember that interactions between particles (or 
holes) and QCSs are suppressed at large angles. Then each patch can be treated incoherently, 
except at boundaries where several patches meet. Knowing the magnitude of the gap allowed 
us to estimate the size of each patch. 

We derived the effective Lagrangian for Goldstone modes by regarding the system as one 
dimensional chains with weak transverse tunnelings. After inclusion of transverse hopping 
effects, the spectrum of Goldstone bosons propagate in three dimensions. As a consequence, 
quasi-long range order in 1 + 1 dimensions becomes, in 3 + 1 dimensions, true long range 
order. Corrections can be systematically estimated through an expansion in Aqcb/pf- 
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Our arguments imply that as the quark density increases, the Fermi surface acquires 
more patches, which have an increasingly large symmetry. Accordingly, there is a series of 
phase transitions as the symmetry of the Fermi surface increases. At the same time, the 
Goldstone bosons become more like those in 1 + 1 dimensions. Then long range order in 
3 + 1 dimensions becomes closer to quasi-long range order in 1 + 1 dimensions. This suggests 
that both chiral and rotational symmetries are effectively restored, but in a manner rather 
different from common assumption. 
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